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ABSTRACT 

We compute the properties of a geometrically thin, steady accretion disk surrounding a central 
rotating, magnetized star. The magnetosphere is assumed to entrain the disk over a wide range of 
radii. The model is simplified in that we adopt two (alternate) ad hoc, but plausible, expressions for 
the azimuthal component of the magnetic field as a function of radial distance. We find a solution for 
the angular velocity profile tending to corotation close to the central star, and smoothly matching a 
Keplerian curve at a radius where the viscous stress vanishes. The value of this "transition" radius is 
nearly the same for both of our adopted £>-field models. We then solve analytically for the torques on 
the central star and for the disk luminosity due to gravity and magnetic torques. When expressed in a 
dimcnsionless form, the resulting quantities depend on one parameter alone, the ratio of the transition 
radius to the corotation radius. For rapid rotators, the accretion disk may be powered mostly by 
spin-down of the central star. These results are independent of the viscosity prescription in the disk. 
We also solve for the disk structure for the special case of an optically thick alpha disk. Our results are 
applicable to a range of astrophysical systems including accreting neutron stars, intermediate polar 
cataclysmic variables, and T Tauri systems. 

Subject headings: accretion disks — magnetic fields — neutron stars — pulsars: X-ray — X-rays: 
binaries — stars: individual ( T- Tauri, FU Orionis) 



1. INTRODUCTION 

In a wide class of objects, from proto-stars to neutron stars, the central accreting object is expected to be endowed 
with a magnetic field of dipolc strength sufficient to influence the motion of matter in the inner parts of the accretion 
disk. The degree of this influence depends on the (unknown) details of the interaction of the dipole field with the 
accreting fluid. It is not clear a priori whether the magnetosphere penetrates the disk, or whether it is capable of 
transmitting significant torques if it docs. In the most commonly accepted model, the stellar magnetic field entrains 
the inner accretion disk as a result of Rayleigh- Taylor and other instabilities and exerts a torque on the disk, whose 
sign depends on the relative angular velocity of the disk and the star. This can result in a spin-up or spin-down torque 
on the star, depending on the value of the inner radius of the disk and on the accretion rate. 

In this work, following the formulations of Wang (1987, 1995), Livio & Pringle (1992), and others, we adopt a simple 
model of distributed magnetic torques on the disk, in which their magnitude depends on the ratio of the local orbital 
angular velocity to the stellar rotation rate. The underlying assumption is that the external magnetic dipole penetrates 
the (thin) accretion disk for a wide range of radii (Livio & Pringle 1992, Wang 1996). The adopted model allows us to 
compute the spin-up/spin-down torques on the central star, as well as the luminosity of the disk, with no reference to 
the viscosity law or to the actual value of the pressure and other disk variables. If, additionally, one adopts a specific 
prescription for the viscosity, e.g., as in an alpha disk, it is also possible to compute the detailed radial profile of the 
thermodynamic variables of the disk. 

The inner termination radius of an accretion disk is of considerable astrophysical interest, as it affects the flow of 
energy and angular momentum in the accretion process. In the early literature, the termination radius of the disk was 
thought to increase smoothly across the corotation radius as the mass accretion rate drops, with matter being ejected 
from the system by the super-Keplerian magnetosphere as soon as the disk is pushed out beyond the corotation radius 
(e.g., Davidson & Ostrikcr 1973, Illarionov & Sunayev 1975). More recently, it has been argued that the termination 
radius does not necessarily become larger than the corotation radius for the so-called "fast pulsars" (Wang 1987), but 
that accretion may be stopped when "the stellar magnetic field imparts more angular momentum to the disk plasma 
than is removed by internal viscous forces in an unperturbed Keplerian disk" (Wang 1995). Other authors agree that 
the disk terminates within the corotation radius, but argue that accretion will proceed even for "fast pulsars," because 
no matter is actually ejected from the disk by the rotating magnetosphere 3 (Spruit & Taam 1993; Rappaport, Fregeau 
& Spruit 2004). Our results support the latter viewpoint. 

The plan of the paper is straightforward. After briefly reviewing some of the more relevant observations, we write 
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3 The actual configuration of the magnetic field can only be determined by time-dependent MHD simulations (e.g., Kato et al. 2001). 
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provide a definitive answer as to what critical "fastness" parameter of the system (see §4, eq. [8]) is required for the onset of the 'strong' 
propeller effect (e.g., Ustyugova et al. 2006). 
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down the governing angular momentum equation, we solve it, and then we discuss the implications of this solution, 
particularly for the torque on the central star and for the disk luminosity. In particular, in §3 we present an overview 
of our magnetically torqued disk solution, in §4.1 we discuss the adopted torque model, and in §4.2 we determine 
the inner radius of the viscous disk. In §5 we compute the torques acting on the disk as well as the concomitant 
luminosity. Our general results are summarized in graphical form in §6, while in §7 we give a detailed solution for 
a magnetically torqued alpha disk. The results are summarized and discussed in §8. An alternate B-field model is 
discussed in Appendix B, previous work is reviewed in Appendix A, and in Appendix C a model of a magnetically 
dominated (zero-viscosity) Keplerian disk is presented. 

2. ESTIMATES AND OBSERVATIONS OF TORQUES. 

Several classes of astronomical sources involve accretion onto a rotating central star supporting a strong magnetic 
field. These include young (proto)stellar objects, such as T Tauri stars, as well as accreting stellar remnants, i.e., white 
dwarfs (cataclysmic variables and polars) or neutron stars (accretion-powered X-ray pulsars). 

It is now possible to directly measure the magnetic field strength in the inner parts of certain proto-stellar accretion 
disks, such as the one around FU Orionis (Donati et al. 2005). Accretion torques have been studied early on for white 
dwarfs (e.g., Lamb & Melia 1987). However, the most detailed information on accretion torques comes from studies of 
X-ray pulsars (Bildsten et al. 1997, and references therein). In particular, it has been found that abrupt transitions 
occur from spin-up to spin-down of the neutron star, with no clear change of the pulsed luminosity (which is presumed 
to originate in matter that is channelled to the magnetic poles), and so perhaps without a large change in the mass 
accretion rate. Nelson et al. (1997) take the view that in Roche-lobe overflow systems the disk may change the sense 
of its rotation from prograde to retrograde, but Li & Wickramasinghc (1998) deem this an unlikely possibility. While 
the transitions remain unexplained, we adopt the conventional view that the disk is prograde, but the magnetosphcrc 
can transmit angular momentum of cither sign between the star and the disk, depending on the location of the inner 
edge of the disk and the mass accretion rate. 

Neither the physics of angular momentum transport by the stellar magnetosphere, nor the mechanism of penetration 
of the magnetic field into the accretion disk is well understood. The main uncertainty concerns the radial extent of the 
disk region threaded by the external magnetic field, as well as the degree to which the magnetic field threading the disk 
is reduced in magnitude by screening currents, and yet there is a remarkable agreement as to the value of the radial 
distance from the star at which magnetic stresses balance hydrodynamic stresses. Various authors, using very different 
physical assumptions, obtain rather similar values for this radius, at least for the case when it is within the corotation 
radius. As we show in Appendix A, because of the rapid, 1/r 6 , variation with distance of the squared dipolc field 
strength, and because of a universal scaling of physical quantities in a thin rotation-supported disk, detailed estimates 
agree, within a factor of < 3, with the simplest dimensionally correct formula 

r m = (GM)- 1/7 M" 2/ V /7 (1) 

(e.g., Lamb, Pcthick, & Pines 1973; Rappaport & Joss 1976; Ghosh & Lamb 1979; Wang 1987; Arons 1993). The 
modest uncertainty in the estimate of the magnetospheric radius is further reduced by a square root when one computes 
the "material" torque on the star 

N m ~ My/GMr m oc M 6/r . (2) 
Observations of X-ray pulsars, indicate a steeper dependence of the spin-up torque on luminosity, 

N~ L 13 (3) 

with (3 > 6/7 (Bildsten et al. 1997), and possibly even (3 > 1 (Parmar et al. 1989). This suggests that either the 
observed luminosity is not proportional to M, or the torque is given by an expression different from eq. ([2]), or both. 

In this paper, we investigate in detail the non-trivial variation with M of the (magnetic and material) torques exerted 
on the star and of the disk luminosity that follow from the adopted model of the magnetosphere-disk interaction. We 
find that at high accretion rate N oc M 9 / 10 (for the alternate model considered in Appendix B, TV oc A'/ 6 / 7 ). At low 
accretion rates, the torques on the central star naturally reverse sign and yield spin-down. While the exact value of 
the spin-up and spin-down torques somewhat depends on the choice of the magnetic torque model, we find it to be 
generally true that at low accretion rates the luminosity of the disk is not proportional to the mass accretion rate. In 
fact, the disk luminosity is greatly enhanced by the torques on the disk whenever the star is being spun down. 

3. OVERVIEW OF DISK SOLUTION 

In the literature, one can find two approaches to the problem of an accretion disk interacting with the magnetosphere 
of a central rotating magnetic dipolc. In the first, a standard thin disk adjusts the rotation rate of its inner edge to 
that of the star — and may be entrained by the external magnetic field — within a narrow boundary layer (e.g., Ichimaru 
1978; Scharlcmann 1978; Arons 1993). In the second, the magnetic field entrains the disk over a wider region (e.g., 
Ghosh & Lamb 1979; Wang 1987; Li & Wickramasinghe 1998; Rappaport, Fregau & Spruit 2004), allowing a smoother 
variation of the torque on the central star with the mass accretion rate. Matt & Pudritz (2005) suggest that one or the 
other approach may be appropriate, depending on the accretion state of the system. For the purposes of this paper, 
we adopt the second approach. 
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The actual profile of angular velocity in the entrainment region is not known, but it is clear that the angular velocity 
has a maximum, at a radius r\. Most authors assume that the disk is nearly Keplerian up to r±, and that radius is 
taken to be close to the magnetospheric radius, T\ ~ r m . The viscous torque is usually taken to vanish only at r±. 
It is understood that this procedure is not consistent, i.e., the angular velocity cannot be both Keplerian (fi oc r -3 / 2 
for r > r\), and have a maximum (dCl/dr = 0) at r\. The inferred magnetic stress at T\ varies by a factor of four, 
depending on whether the former or the latter is taken to hold (e.g., Li & Wang 1996), because the magnetic stress is 
taken to be proportional to d(£lr 2 )/dr (see eg. [To]) . 

In this paper we adopt a model for the magnetic torques on the disk and then proceed to solve the azimuthal 
component of the equation of motion. Wc assume that the viscous torque vanishes at a certain radius ro, and that 
the disk is Keplerian at least for all r > tq. This implies that ro > r%. The value of ro follows from the adopted 
model of the magnetic torque. The viscous torque can vanish at ro because the magnetospheric stresses are sufficiently 
high to remove angular momentum (from a thin annulus at this radius) at the rate required to sustain the prescribed, 
constant and uniform mass accretion rate. The magnetospheric stresses on the disk can only increase inwards from r . 
Accordingly we take the viscous torque to vanish for all r < r as well. No other assumptions are needed to solve for 
O(r), once the magnetic torques are specified. It is easy to find an angular velocity profile that smoothly matches the 
stellar rotation rate at small radii and the Keplerian curve at ro f ^4.2[) . 

An additional point should be mentioned. If the stellar magnetic field entrains the disk, it exerts a torque on it, 
i.e., it either deposits angular momentum in the disk, or removes some angular momentum from it. This affects the 
amount of angular momentum transported in the disk by the viscous torques, and hence the rate of energy dissipation 
and the luminosity of the disk ( £!5.2p . 



4. MAGNETICALLY TORQUED THIN DISKS 

4.1. Angular momentum transport 

We start by defining a few fiducial radii within the accretion disk which we will utilize in this work. The corotation 
radius is taken to be the radial distance at which the Keplerian angular frequency, Qk, is equal to the rotation 
frequency of the central star, uj s : 

1/3 

(4) 



= ( GM \ 



where M is the mass of the central star. We also define a convenient dimensional magnetospheric radius: 

r m = (GM)- 1/7 M- 2/ V /7 , (5) 

where M > is the steady-state mass accretion rate, and /i is the magnetic dipolc moment of the central star. Here 
and throughout this work we take eq. © to be a formal definition of r m , however, our value of r rn is identical to the 
inner disk radius in the model of Arons (1993). Because the ratio of these two fiducial radii appears quite often, we 
define a dimcnsionless parameter, £, which relates all four system parameters of the problem (M, M, [i, u> s ): 

^— ■ (6) 

r c 

For a fixed system, this parameter is a measure of the mass accretion rate: £~ 7 / 2 oc M , for M = const, fi = 
const, uj s = const , so we may define a fiducial accretion rate Mq through 

£ -7 / 2 = M/Mq , for M — const, // = const, u> s = const . (7) 

Finally, we define ro as the radial inner boundary of the region where the accretion disk is Keplerian and the viscous 
torque is non-vanishing. The "fastness" parameter corresponding to this transition radius is defined as 



^if(r ) 



, 3/2 



(8) 



The equation governing conservation of angular momentum for a thin accretion disk subjected to distributed magnetic 
torques from the central star is: 

* (ttr 2 )=r visc + r B <o , (9) 



AuHr dr 

where r is the radial coordinate, f2 is the local orbital angular frequency of the disk material, 2H is the full disk 
thickness (which is a function of r), r v i SC is the viscous torque per unit volume, and Tb is the magnetic torque on the 
disk per unit volume. Here, the advection of angular momentum by matter flowing inwards is driven by the viscous 
and magnetic torques, where a positive sign for Tb indicates that angular momentum is deposited in the disk. The 
vertically averaged viscous torque per unit volume can be written as: 

T ^ = ^ri^ 2H ^ = ^rlM > (10) 
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where t rt j, is the vertically averaged- (and T rt j, is the height integrated-) r - <j) component of the viscous stress-energy 
tensor; the two are related through T r( p — t r< f, ■ 2H. 

We take the magnetic torque on the disk to be distributed over a range of radial distances with the magnetic torque 
per unit volume given by: 

_ r d{B zB4> ) 

rB "47^^' 
so that the vertically averaged magnetic torque is 

±[r Bdz = ^, (12) 

2H J 4ttH ' V ' 

where, in the last term, the magnetic field is evaluated at the top "surface" of the disk, i.e., at z = H. 

Without a full, self-consistent magnetohydrodynamic solution of the disk equations, there is no simple way of 
determining B. In principle, the magnetic torque on the accretion disk should be computed self consistently, but for 
purposes of gaining insight into the global effects of magnetic torques we prefer to adopt a reasonable analytic model 
rather than working with a more complex problem that is not necessarily any more valid. To obtain results that 
would be applicable to any thin disk model, we simply adopt the following somewhat ad hoc but physically plausible 
prescription for at z = H (Livio & Pringle 1992; Wang 1995; Rappaport, Fregeau, & Spruit 2004): 

B4, = Bl~B,\l-?£j , (13) 

which we take to hold for r > r c (see also, Wang 1987; 1995; 1996 for more sophisticated and physically motivated 
versions of this relation). Here and in the following we assume, without loss of generality, that B z > 0. The expression 
for B^ given in eq. (|13jl has the property that the magnetic torque vanishes at the corotation radius, while B^ is 
comparable to B z at large distances. The latter is compatible with the fact that B^ cannot be larger than B z 
over significant distances, for reasons of equilibrium and stability of the field above the disk plane (see discussion in 
Rappaport et al. 2004). As the field is wound up by differential rotation between the star and the disk, the azimuthal 
component first increases, but the increasing energy in the azimuthal field component pushes the field configuration 
outward into an open configuration when the azimuthal component becomes comparable to the poloidal component. 
This was proven by Aly (1984, 1985) in a rather general context, and worked out in some detail for the case of a disk 
around a magnetic star by Lynden-Bell & Boily (1994). For radial distances inside the corotation radius (r < r c ) we 
adopt one of two plausible expressions for B§: 



B{~B z {l- — ) or Bi 1 - 



where the first version, B^, is just an analytic continuation of the expression given in eq. fj 13|) . while the second 
version, B^ 1 , has the aesthetic advantage that \B,p\ — > \B Z \ as r — * 0, as opposed to becoming progressively S> B z (for 
Keplerian Vl). However, B^ 1 has the disadvantage that the functional form for the magnetic torque must be switched 
at r = r c . We utilize both of these functions in this work, but give results in the text for B^ only since these are 
algebraically less messy. The corresponding results for B 1 ^ are given in Appendix B. We note that a form similar to 
B 1 ^ was discussed extensively by Erkut & Alpar (2004), who also derive an equation very similar to our eq. (|B5|) . 

These two prescriptions (cq. [14] ) for the azimuthal field can be discussed in terms of the magnetic diffusion coefficient 
77. As argued by Livio & Pringle (1992) the azimuthal field is created by vertical shearing. The poloidal field is stretched 
by a factor of 2-KrjB in a time 2tt/\H — uj s \, i.e., the azimuthal field is created at the rate (r/H)\B z (Q, — u> s )\. In steady 
state this has to be balanced by the diffusion rate rjcPB^/cPz ~ riB^/H 2 , resulting in B^ ~ rHB z {uj s — Vl)/r). Campbell 
(1992) derived the last relation from the induction equation, while Wang (1987) obtained it from considerations of 
buoyancy. Our two prescriptions for B§ then correspond to r\ ~ rHuj s , or r\ ~ rHfl, respectively. Close to the 
corotation radius both prescriptions coincide. In the second prescription, the diffusivity is larger by a factor r/H than 
the largest possible value of the coefficient of kinematic viscosity in subsonic turbulence, v max ~ c s H ~ H 2 fl ~ r]H/r, 
where c s is the speed of sound, i.e., we do not necessarily assume that the magnetic diffusivity and the disk viscosity 
are of the same origin (cf. Campbell 2000). 

We note that Agapitou & Papaloizou (2000) use a diffusivity prescription resulting in a torque that differs from 
our prescription by a constant factor C (their eq. [14]). Their discussion suggests that a large value of C is more 
physical than the value C = 1. However, the effect of smaller assumed diffusivity seems to be offset by the effects of 
an inflated field (see below). We further note that in the following, a rescaling of the diffusivity by a constant factor 

l/7d corresponds to rescaling of rj/ 2 by jd, and of \i by ^/j d - Wang (1995) showed that the torque on the central star 

is independent of 74, when expressed in units of My/ GMro (see cq. [37], where we have neglected the contribution 
-lo s R 2 s in eq. [55]). 

With the azimuthal field given by B^ in eqs. p^]) . (fT4"|) the complete height-averaged angular momentum equation 
now reads: 

M d ln s Id,— n\ B 2 r 

(Mr 2 ) = (Trd>r 2 ) H — (1 - 9, uj s ) . 

4irHrdr K > 2Hr dr V ' AttH V 1 s ' 

(15) 
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Fig. 1. — Ratios of ro, the inner radius of the viscous disk, to the 
corotation radius, r c , and to the magnetic fiducial radius, 7*771 , as 
functions of the parameter £ _1 = r c /r m . The results are shown 
for the two different prescriptions for (see text). The vertical 
line indicates the v alue of l/£ corresponding to zero torque on 
the central star (eq. 1360 . The horizontal dot-dashed curve is the 
Ghosh Bz Lamb (1979) value for ro/r m . 
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Fig. 2. — The matching solutions in t erms of Q/lu s in the sub- 
Keplerian regime as found from eq. 1161 . The curves are labeled 
by the values of the parameter £. 



We consider a central dipole aligned with the stellar spin axis, and an accretion disk that is axisymmetric and 
perpendicular to the spin axis. In the following we will assume that the poloidal field is given by the dipole formula 
B z = [i/r 3 . We expect that the results obtained in this paper will not differ qualitatively if this assumption is relaxed. 
For example, the expression of Arons (1993) for the inner radius of the disk, r m , which we reproduce in eq. (|A12p of 
Appendix A, is based on an exact solution of Aly (1980) for the field structure in the presence of a diamagnetic disk 
(in which, close to the inner edge of the disk the field is enhanced by a factor of r/H), and yet, it agrees with our 
result for ro: ro/r m ~ 0(1) for r m < r c ; see eqs. (|17|) . (|B7[) . and Fig. [1] Likewise, considering numerical solutions of 
realistic disks, Bardou & Heyvaerts (1996) and Agapitou & Papaloizou (2000) show that the field lines may be inflated 
outside the corotation radius, leading to a reduction in the value of the poloidal field relative to the dipole value that is 
assumed here. However, this does not seem to strongly affect the results discussed in this paper. Eventually, Agapitou 
& Papaloizou (2000) find values of torque which are the same as ours, within a factor of a few [see our H5.3] , 

4.2. Matching solutions in the sub-Keplerian regime 

Unlike most other workers, to describe the steady disk we use the same disk equation for all radii (eq. [15]). Since 
the model torques depend on the radial profile of the angular velocity in the disk, and the angular velocity profile 
depends on the torques, the azimuthal equation of motion of the accreting fluid must be solved in a self-consistent 
manner. Assuming no mass loss from the disk, we find a solution which asymptotically matches the rotation rate of 
the star at small radii, and smoothly matches a viscous Keplerian disk at a certain radius, ro. The value of this radius 
also follows from the same eq. (fTS"|) . 

We begin by arguing that the disk may be taken to be Keplerian where the viscous torque is non-vanishing, with 
the corollary that T r< f, = for r < ro. In a steady, thin accretion disk, the radial pressure gradients are a factor of 
(H/r) 2 smaller than gravity (Shakura and Sunyaev 1973; H/r << 1 is the dimensionless disk thickness), and hence, 
in the absence of strong external disturbance the disk is rotation supported, i.e., nearly Keplerian everywhere. On the 
other hand, if the disk extends within the corotation radius, its angular velocity must reach a maximum (at a smaller 
value r±) before it can match the lower angular velocity of the star. At the maximum, the viscous torque must vanish. 
We expect the viscous torques to be zero also for all r < r\\ this is consistent with the fact that magneto-rotational 
instability (Balbus & Hawley 1991), which is thought to be responsible for the presence of an effective viscosity in the 
disk, does not operate in the region of radially increasing angular velocity dfl/dr > 0. Hence, the viscous torque 
vanishes for all r < r2, for a certain r% > r\. 

Accretion can proceed in the region r < T2 only because the magnetic torques are sufficiently high to remove angular 
momentum at the requisite rate. In fact, by the argument in the previous paragraph, in the sub-Keplerian region the 
magnetic torque must be dominant, the viscous term in eq. (|15|) only becoming important when \B 2 r (1 — Q,k/oj s ) \ < 
| — MQ,k/2\, where £Ik is the Keplerian orbital frequency. To keep the algebra clean we assume that r-i > ro, i.e., the 
viscous torque vanishes already in the Keplerian region. Since eq. (|15|) does not admit Keplerian solutions for T r( j, = 0, 
we must then have ?*2 = ro- This is not a very restrictive assumption. If, instead of assuming T r( p(ro) = 0, we were 
to take the viscous torque to remove one-half of the angular momentum at the matching boundary, d(T r< pr 2 ) / dr\ ro = 
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— {M / Air)d{Q,Kr 2 ) / dr\ TQ , the estimate of ro in eqs. (fT"7|) - (p~9|) would change by a factor of 2 1 / 5 for £ << 1, i.e., by less 
than 15%, and not at all in the limit £— > oo. In short, we assume T r ^ ^ if and only if r > r$. 

The equation for angular momentum transport (eq. |15| ) in the region r < ro, where viscous stresses are zero, is our 
starting point for exploring formal flow solutions which make the transition from Kcplcrian rotation at ro to corotation 
with the central star as r — > R s : 

M d , 9N B 2 r j . i . , 

^ 2 =--f=(l-f2/u; s ) . (16) 



AnHr dr y ' 4ttH 

This being a first order ordinary differential equation, a single initial condition (at a fixed radius) specifies the solution. 
However, using the freedom of choosing the matching radius, it is possible to match both the slope and the value of 
the solution with a Kcplcrian one. Indeed, we first find the radius ro at which the solution matches Kcplcrian rotation 
by substituting f^r/2 for the derivative on the left hand side, and then use f2(ro) = £Ik(to) = y GMj r^ as the single 
initial condition to solve for f2(r) for all r < r . 
With the described substitution, cq. p6[) reduces to an algebraic equation for ro: 




£7/2^-10/3 (17) 



This is essentially the same as the expression for the inner edge of the disk found by Wang (1995), and also used by 
Kenyon et al. (1996). For any choice of value for the dimensionless parameter £ this equation can be solved numerically 
for the inner radius, ro in units of r c . The asymptotic limits for large and small values of £ are: 

r£ „ 2 l/5£7/10 for f<1 (1 8 ) 

r c 

^~i-±r 7/2 for e»i • (19) 

r c J 

Note that ro < r c , always. Plots of ro/r c and ro/r m are shown in Fig. [T] For a given system, the high accretion 
rate limit of eq. (fTT)) is r oc Af ~ 1 / 5 , by eqs. (J7J), {T5} . Of course, the disk cannot penetrate the surface of the star or 
extend far within the marginally stable orbit (ISCO) predicted by general relativity (e.g., Kluzniak & Wagoner 1995), 
so eqs. (|17M19[) and the solutions discussed below are valid only for ro >> max(i? s , r ms ). 

For any non-vanishing mass accretion rate, M ^ 0, eq. ()17|) allows us to eliminate the unobservable quantity r m 
from all the equations, so that the ratio ro/r c is the only remaining parameter in the expressions for various physical 
quantities. Specifically, the first equality in eq. (fTT)) can be rewritten as 

\ 7/2 

2-] = - ■ (20) 
ro J 2(1 - w) [ ' 

The limit of no magnetic field corresponds to /i 2 = 0, r TO = 0, £ = 0, u> = 0. The limit of M — ► corresponds to 
u> — > 1, and more specifically M/(l — uS) — > 2Mo- 

The counterpart of eq. (fT7|) for our (alternate) prescription (eq. |B7| in Appendix B) has also been derived by 
Matthews et al. (2005). However, their interpretation of ro is different from ours. For Matthews et al., ro (their R t ) 
is the inner edge of the magnetically dominated disk, where the surface density vanishes (i.e., S = for r < R t ), and 
where the radial inflow velocity is large. In our approach, ro is the outer edge of the magnetically dominated inner 
disk; the disk is Keplerian at ro, and therefore the radial velocity must still be much smaller than the value of rO#- at 
this radius. Furthermore, there is no physical reason for the density to vanish at ro. 

The specific angular momentum of the matter flow, I = fir 2 , is governed by the equation valid for all r < ro 
(essentially a rewritten version of eq. |16j): 

Md£ „, /„ t \ 

— T = - B * r \ 1 2 . 21 

r dr \ uj,r z J 



which can be cast in completely dimensionless form by using rescaled variables C = £/\JGAIr c , and 1Z = r/r c 



dK K 4 V ^ 2 . 

This equation admits to analytic solutions which involve incomplete gamma functions, but instead of reproducing the 
formulae we exhibit the solution graphically, after simply integrating this equation numerically. The results are shown 
in Fig. O together with the Keplerian curve. Note that the solutions smoothly match a Keplerian Q(r) at ro, reach a 
maximum at a lower value of n, and finally decrease smoothly until they level off asymptotically toward u> 8 . 

Neither the value of the maximum angular frequency, fi(ri), nor the value of the radius r\ at which this maximum 
is attained, plays any role in the following considerations. This is because the viscous torque is non-zero only for radii 
r > ro > r\. There is no viscous dissipation, and no angular momentum can be transmitted upstream through the 
fluid, for any r < ro. Once the accreting fluid passes through the (imaginary) cylinder at ro, all its angular momentum 
is transmitted to the star, at the rate M^/GMtq. 



Magnetically Torqued Accretion Disks 



5. TORQUES AND LUMINOSITY 

Now that we have found a solution in which the viscous torque vanishes at ro and the disk is Keplerian for all r > ro 
we can use eq. (|15|) to compute the viscous torque in the disk, and the luminosity of the disk. We are also in a position 
to compute the total torque on the star. These results follow from our model assumptions alone, and are independent 
of the detailed structure of the thin accretion disk (equation of state of the fluid, the viscosity prescription, disk height, 
etc.). 



If we integrate both sides of eq. (fl 
of the viscous accretion disk) , we find 



5.1. Viscous Torque 
with respect to r, starting at ro, the largest radius where T r( , 



M (fir 2 - IVo) = 27rT r , 



n 



dr, 



(inner edge 



(23) 



where for the z component of the magnetic field we have taken fi/r 3 to represent an unscreened dipole field whose axis 
is aligned approximately along the spin axis of the central star and perpendicular to the accretion disk. Finally, we 
can carry out the remaining integral and solve for the vertically integrated r — <p component of the stress tensor at an 
arbitrary radial distance: 



T, 



MsfGM 
2nr 2 



18ttH 




(24) 



Here we have made use of our result from Section 14.21 that the disk matter is in Keplerian orbital motion everywhere 
for r >tq. 

5.2. Viscous Heating of the Disk 

We now utilize our result for the vertically integrated stress to compute the local viscous heating in the disk. The 
viscous heating per unit volume is given by Q = tr^rdtt/dr, and the corresponding viscous heating per radial interval 
is dL v - ls /dr = 2irr 2 T r< f 1 dfl/dr. From this we find: 



dL v 



dr 



3GMM 
2r 2 



1- J — 



'VgaI 



Q r 5/2 




(25) 



The first term on the right hand side of the equation is the usual luminosity that results from the redistribution by 
viscous torques of the released gravitational potential energy (Shakura & Sunyaev 1973). The second term represents 
the viscous dissipation of mechanical energy that is input to the disk via the magnetic torques. By utilizing the 

7/2 



definition of eq. ©, r„ 
dLs, 



h 2 /(mVgm) 

3GMM 



we can cast eq. (|25j) into a somewhat more aesthetically pleasing form: 

J/2 



dr 



or 



dL v 



2r 2 



3GMM 




9r l/2 r 3 



3-3- 



dr 



2r 2 



1 + — 



-1 



1 

18(1 -uj) 



3w- 3u-± - 2 




(26) 



(27) 



Finally, eq. (f2"6"| can be analytically integrated to yield the total power released by viscous processes in the disk: 



GMM 
2r 



9 V r 



7/2 




(28) 



The first term represents the gravitational power that would be released in the disk in the absence of magnetic torques. 
Clearly, in this model the disk may be powered in part by the central star — this occurs whenever O(ro) < (4/3) u> s , 
or, in terms of the fastness parameter, when cu > 3/4. In general, after the constraint of eq. ([20]) is used, eq. I|28p may 
be rewritten as 



GRIM (15 - Ulo) _ GMM (15 - 14w) 
~ 2r c 



-2/3 



(29) 



2r 18(1 -u) 2r c 18(1 - u) 

This quantity is positive for all lu < 1. 

The total mechanical power injected into the viscous disk and the part of the magnctosphcrc entraining it at r > ?*o 
is given by 



Ltnt — 



GMM 
2rn 



EL 

9r 3 




(30) 
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Fig. 3. — dL vla /dr vs. the radial distance (in units of r c ) for 20 
values of the parameter £ in equal logarithmic steps over the range 
0.1 < £ < 10 (see eq. 1261 ); with £ increasing from left to right. 
In terms of the fastness parameter, this range of £ corresponds to 
0.113 < w < 0.99984. In these plots dL vls /dr has been normalized 
in terms of GMM/(2r%). 



Fig. 4. — The quantity (dL vla /dr)(r /ro) vs. radial distance (in 
units of r c ) for 20 values of the parameter g i n equal logarithmic 
steps over the range 0.1 < {; < 10 (see eq. 1261 ); with £ increasing 
from left to right. In terms of the fastness parameter, this range of 
§ corresponds to 0.113 < ui < 0.99984. The ordinate in these plots 
has been normalized in terms of GM M / (2ro); the factor r/ro 
renders the plotted functions roughly equivalent to the luminosity 
per natural logarithmic radial interval. 

where the second term on the right hand side is the rotational frequency of the central star times the total magnetic 
torque on the disk (see eq. [53] below). For M ^ 0, this expression reduces to: 



Eq. (f3l"j) can be rewritten as 



GM M 
2r 



GMM 
2r 



1 w(3w - 2) 




9 (1-w) 



GMM 
2r, 



1 w(3w - 2) 



9 (1-w) 



,-2/3 



The difference between the expressions in eq. (|3"2l and eq. (|29p is 



/ GMM \ oj- 2 / 3 (3 - 8lj + 6UJ 2 } 
\ 2r c J 18(1 -u) 



(31) 



(32) 



(33) 



a positive quantity for all u> < 1. 

It is not clear what fraction of the difference L to t — L V i S is released in the disk, and what fraction is released in the 
magnetosphere (perhaps as non-thermal radiation). In the limit of M — * 0, L to t = 2L V i S , and L V is = ^sT c /2, where 
r c = /i 2 /(9r c 3 ). 

5.3. Magnetic torques on the disk and the neutron star 

The total magnetic torque on the accretion disk, rs^isk, is given by the magnetic contribution to —2itr 2 T r ^ evaluated 
as r — > oo (this is equivalent to the integral on the right hand side of eq. |23| integrated to oo). This quantity can be 



(34) 



found from the second term on the right side of eq. (|24|) . multiplied by 2nr , and in the limit r — > oo: 



Tfi.disk 



JJ_ 

9r 3 



3-2. 



For M ^ eq. (f34|) may be rewritten in terms of the fastness (eq. [8]): 

rs.disk = [m^GMV^ • (35) 

It may be interesting to note that the net torque on the disk vanishes, i.e., Tg^isk = 0, when r2(r ) = (3/2)cj s . In the 
limit of M 0, rs.disk -» t c = ^ 2 /(9r 3 c ) = M VGMr~ c /9. 
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Fig. 6. — The integral curves for dL vla /dr as given by Fig. 3. 



Fig. 5. — The quantity dL v - m /dr vs. the radial distance (in 
units of r c ) for 20 values of the parameter g i n equal logarithmic 
steps over the range 0.1 < § < 10 (see eq. 1261 ); with £ increasing 
from left to right. In terms of the fastness parameter, this range 
of £ corresponds to 0.113 < lu < 0.99984. These plots have been 
normalized to dL vlB /dr for a Shakura-Sunyacv disk with the same 
values of M and inner radius, tq. 

Here, we can compare our results with the numerical work of Agapitou & Papaloizou (2000). If the disk extended 
only between ro = r c and r = r ox t, as Agapitou & Papaloizou (2000) assume in their calculations, we would have 
obtained 



n 



(/x 2 /3r c 3 ) 1/3 - rl/r 3 cxt + (2/3)(r c /r cxt ) 9 / 2 



instead of eq. ([34]) . Agapitou & Papaloizou take = CB^ (see our eq. [13] ) with B z found numerically to be reduced 
relative to the dipole value that we assume. The torque N(C) they find in their numerical solution is a few times 
larger than n, with N{C)/t\ ~ 2 for C ~ 100 — see eq. [26] and Fig. 14 of Agapitou & Papaloizou (2000). Thus, our 
simple model yields results for the torque which are not too different from those found for more realistic magnetic field 
configurations. Since in reality ro ^ r c , we note that rescaling of the product of the toroidal and poloidal magnetic 
fields, B^Bz, by a constant factor, Ci, affects our results only to the extent that ro varies, so that the values of ui are 
rescaled (eqs. [T7], [20] )■ At the same time the constraint of eq. (fTT]) guarantees that there is no explicit dependence 
on C\ in eq. (|35[) . and in many of our other equations. Formally, ro and TB,disk arc invariant under the rescaling 
Bj,B z C 1 B 4> B Z and M CfM, see eq. ©. 

The magnetic torque on the central star is just the negative of eq. (j34|) . To find the total torque on the central star 
we add the material torque: 




T CS = M (jGMr Q - co s i? 2 ) - £3 

(36) 

where R s is the radius of the central star. (Note that uj s R 2 s is taken to be negligible in the production of our graphical 
results.) For high accretion rates, the material torque is the leading term in eq. (f3"6"]l . Teg ~ M^/GMtq oc M 9 / 10 . 

If we combine eqs. (|36p . ([S]), and (|20p , we can find an analytic expression for the condition that there is zero torque 
on the central star (under the approximation that we neglect the oj s R 2 term in eq. [36], cf., eq. |39|). We find 
r /r c = (20/21) 2 / 3 , and r m /r = 10 2 / 7 , yielding £ T=0 = 1.869 and r r=0 = 0.9680r c = 0.5179r^=°. The corresponding 
values for our alternate model (Appendix [B| are very similar: £ r= o = 1-869 and r T=0 = 0.9666r c = 0.5172r^„ =0 . 
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TABLE 1 
Some fiducial radii 



eq. (8) 




eq. JT71 




cq. 






ro/r c 




r 1 


M/M 


comment 


0.9524 


0.9680 


1.869 


0.5351 


0.1120 


res = 


0.7808 


0.8478 


1 


1 


1 


definition of Mq 


3/4 


0.8255 


0.9269 


1.079 


1.304 


L viB = GMM/(2r ) 


2/3 


0.7631 


0.7634 


1.310 


2.576 


r B = 0, Ltot = GMM/(2r ) 



6. RESULTS IN GRAPHICAL FORM 

The inner disk transition radii, r , as given by eqs. (fT7|) . (|B7[) . are shown in Fig. [1] The four plotted curves are ro 
in units of r m , and ro in units of r c , for both of our prescriptions for B^. These radii are plotted as a function of the 
parameter = r c /r m . Some fiducial dimensionlcss values of ro and r m can be found in Table [TJ 

In the slow rotator limit (i.e., weak B or high M; large £ _1 ), ro/r c — > 0, as expected, while ro/r m asymptotically 
approaches 2 2 / 7 = 1.22 for the B¥ model, and ro/r m is larger for the B^ model, slowly increasing with increasing 

In the opposite limit of a rapid rotator (i.e., strong B or low M; small we see that ro — > r c as £ _1 — > 0, as 

discussed extensively by Rappaport ct al. (2004). In that paper, the authors suggest that for reasonably fast central 
rotators, the matter is not ejected via a propeller mechanism, but rather the accretion disk adjusts its structure to 
allow matter to reach r c (see also §4.2 in this article). Note that for this and other reasons, ro/r m varies quite strongly 
with r c /r m , in contrast with the Ghosh & Lamb (1979) model, where the same ratio, ro/r m , is a constant equal to 
about 0.4. In particular, we find that ro/r TO — * as — ► 0. 

In Figures 3 through 5 wc plot the viscous luminosity in the disk per radial interval, dL wls /dr, eq. (|26p . with three 
different normalizations. In the first of these (Fig. 3), dL V i S /dr is normalized to the differential gravitational luminosity 
at ro, viz, GMM /ftr 3 ). Fig. 4 is dL vla /dr multiplied by a factor of r/ro to yield the viscous luminosity per natural 
logarithmic radial interval. In Fig. 5, dL vls / dr is normalized to the same quantity for a Shakura & Sunyaev disk. In all 
cases, the pattern of the curves indicates how the inner radius of the disk moves outward as the parameter £ increases 
(i.e., as the central object becomes a faster rotator), and ultimately asymptotically approaches r c . For very rapid 
rotators, the plots clearly indicate that the disk energetics are dominated by the mechanical energy being pumped in 
by the magnetic field, and dissipated (non-locally) by viscous stresses. For the rapid rotators, the extra energy being 
pumped into the disk is powered by a spin-down of the central star. Finally, in Fig. 6 wc show the integral viscous 
disk luminosity profiles, i.e., L vm released for radii > r. 

We note that in Figs. 3-6, due to the way the curves have been normalized, the large disk luminosities for rapid 
rotators are only with respect to L v ; s in the absence of magnetic torques. The absolute values would steadily decrease 
as £ grows, since M is decreasing in the normalization factor GMM /(2ro). The absolute value of the disk luminosity 
(i.e., not normalized to 1/ro) is shown in Fig. 7 as a function of M . In order to make these plots as generic as possible 
we have normalized the luminosities to GMMo/{2r c ), where r c is fixed for a given system, and M is normalized to 

i.e., a value that is also fixed for a given system. The three luminosities shown in the figure 
are (i) L v is, as given by eq. (|28p . and which includes the mechanical energies put into the disk via the gravitational 
field and the magnetic torques; (ii) = GMM /(2r ) which results from the deposited gravitational energy alone; 
and (iii) the total energy released in the system, including that due to an equivalent "frictional" dissipation via the 
magnetic field (see cq. [31]). The effects of magnetic heating of the disk are very apparent for low values of M, i.e., 
the fast rotator case. The effects of magnetic "cooling" , for slow rotators, are much less apparent. 

The torques on the central star are shown in Fig. 8. The four curves are combinations of |tb| and | Ttot | vs - -^vis 
or Ltot- Here L vls and L to t are defined by eqs. (|28|) and ([3Tj) . respectively. The torques, tb and r to t, are defined in 
cqs. (|34]l and (|36|) . respectively, where r to t = tcs + Muj s R 2 s . For high values of the disk luminosity we see that the 
magnetic torques are typically an order of magnitude smaller than the matter accretion torques. At lower luminosities, 
the magnetic torque switches sign at a luminosity of about (1.2 - 1.3) GMM / (r c ), while the total torque (including 
that due to accretion) reverses sign only at luminosities about an order of magnitude lower (cf. Table [T]). Note that 
the negative torques tend to a constant value, as the disk luminosity decreases (i.e., as the mass accretion rate drops). 

In terms of the spin- up/spin-down dependence on the mass accretion rate, our results are qualitatively similar to 
those of Wang 1987, 1996, and quite different from those of Ghosh & Lamb (1979). With the help of cq. ([IT]), and 
neglecting the term uj s R 2 , the torque of eq. (|36|) can be expressed as a dimensionlcss function gi of the fastness 
parameter: 

r cs = (Kly/GMro) gi(w), (37) 

or, with a different normalization, 



T CS = (AIy/GMr c ) uj 1/3 gi (u), (38) 
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Fig. 7. — Luminosity as a function of M. The reference mass 
transfer rate, Mo, is the value of M which, for given values of 
fi, M, and r c , yields a ^-parameter value of unity. The curve 
labeled L VIa desc ribe s the total viscous luminosity of the disk, 
as given by eq. 1281 . while the curve labeled Ltot is the total 
mechanical energy put into the system via the magnetic field 
and grav ity - including the magnetic "frictional" dissipation term 
(eq. 1311 ). The curve labeled L^ f is the gravitational energy re- 
lease, GMM/(2ro). All luminosities are expressed in units of 
GM Mo/(2r c ). In this paper, we do not discuss the contribution 
to luminosity of the power released by matter in transit from the 
viscous disk to the stellar surface. 



Fig. 8. — The torque on the central star due to magnetic 
torques (right set of curves, ne gati ve of eq. 1341 ) and mag netic 
plus matter torques (left set, eq. |36|). vs. L vla and Ltot (eqs. |28j . 
[31|). The M dependence of these luminosities is shown in Fig. 7. 
The torques and luminosities are normalized to Mq \J GMr c and 
GMMo/(2r c ), respectively, where Mq is defined in the caption to 
Fig. 7. Only the absolute magnitudes of the torques are plotted 
to allow for the use of a log scale. 



with 



10 
~9 



1 



1 



20 1 -u> 



(39) 



As in Wang's (1996) dimcnsionlcss torque, gj is a fairly flat function up to the zero at u = 20/21 = 0.9524 « 
(0.9680) 3 / 2 , in general agreement with the observations of Finger et al. (1996) and in contrast to Ghosh & Lamb's 
(1979) steeply decreasing torque which vanishes at u « 0.35. Remarkably, the corresponding function gu in our 
alternate model, although also fairly flat with a value higher by a factor ~ 63/60, has a first zero at a very similar 
value of w = 0.9502 (eq. [B15] ). 

As discussed in § 15.31 above, the condition for zero torque on the central star is £ T= o = 1-869 (£~i = 0.535) for 
both of our models. This corresponds to ro/r m ~ 0.52, a value which coincidcntally happens to be in approximate 
agreement with the Ghosh & Lamb (1979) model prediction of ro/r m ~ 0.42 for any torque (Fig. [1]). However, the 
value of fastness at this point, uj t= q = 0.95 is quite different from the Ghosh & Lamb (1979) value of u> « 0.35. 

The zero torque constraint on £ leads to a relation between the equilibrium spin of a neutron star at fixed M and 
the surface B field: P s ~ 2.9£ t Jq 2 fj^Y 'M 18 3 ^ 7 ms, where ^27 is the neutron star magnetic moment in units of 10 27 G 



cm , and Mis is the mass accretion rate in units of 10 g s , or approximately the Eddington limit for a neutron 
star. The equilibrium spin period turns out to be 

P s = 1.14 ms (M 18 )- 3 / 7 B 6 g /7 , (40) 

where Bg is the surface value of the magnetic field in units of 10 9 G. Unfortunately, at present it is not clear how close 
to equilibrium the various X-ray pulsars are. 

Observations of X-ray pulsars indicate a bimodal behavior of spin-up and spin-down, with a rapid change between 
two states of similar |P| magnitude, but of opposite sign of the rate of change of the pulsar spin, P (Bildsten et 
al. 1997). The change between the two states can be quite abrupt, and occurs without a large change in (pulsed) 
luminosity. In our model, a change in luminosity over a fairly narrow range of luminosity could, indeed, lead to a 
change from spin-up torque to a spin-down torque of comparable magnitude (or vice versa), see Fig. 8. Further, the 
spin-down torque quickly tends to a constant value, as the luminosity drops. This is in rough qualitative agreement 
with the observations, however the model seems to admit a wider range of torques, particularly of spin-up torques, 
than is oberved. In particular, recent observations suggest very little change, if any, of the magnitude of the torque 
during the observed transition from spin-up to spin-down in 4U 1626-67 (Krauss et al. 2007, see their Fig. 5). 
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7. COMPLETE THIN ALPHA-DISK SOLUTIONS 

Our expression for the vertically integrated viscous stress (eq. [53]) can be used to derive the full set of thin accretion 
disk variables if we specify a prescription for the stress. The most often utilized of these is the "alpha-prescription" 
of the Shakura-Sunyaev model (1973) in which \T r A is taken to be aP ■ 2H, where P is the disk pressure in the 
midplane, and a is a dimensionless parameter which specifies the strength of the viscous forces. Equation (|24[) then 
directly yields: 

2P ff= ^xF( v / ^, w ) , (41) 

with 

Here, the first term in parentheses in the expression for F is the usual factor /(r) 4 in the Shakura-Sunyaev solution. 
If it were not for the constraint of eq. (|17[) . putting r m = one would recover the Shakura-Sunyaev solution, but 
this would be inconsistent as r is no longer an arbitrary constant, and instead must tend to zero with r m . After 
eliminating r m from eq. (|4"2"]) with the help of eq. (|2T)]) . we obtain 

F(y, uj) = [18(1 - lj)]- 1 (18 - 20y - 18a; + 2luoy - 3uy 7 + 2y w ) , (43) 

in the domain 0<y<l,0<o;<l. Here, y = yro/r, and u> is defined in eq. ([8]). The limit of vanishing field, p, — > 
is given by 

F(y,0) = (l-y)-(y-y 10 )/9, (44) 

because u> = in this limit (eq. |17j). For a fixed value of r (or one bounded from below), y — > 0, F — > 1, in the same 
limit. 

Three of the four remaining SS-disk equations remain unchanged. For the equation of state we take: P = pkT/m, 
where m is the mean molecular weight per particle. The vertical force balance equation is given by: P ~ GM pH 2 r~ 3 . 
Radiative transport in the vertical direction is represented by: X" 4 ~ npHT^, where T and T e are the midplane and 
effective temperatures of the disk, respectively, and k is the radiative opacity evaluated at the disk midplane. Finally, 
we retain the fifth SS-disk equation for the heat dissipation per unit surface area of the disk: T — QH ~ aHPfl = <tT^. 
In doing so, we have explicitly neglected the magnetic "frictional" heating; this latter heat input is never dominant 
over the mechanical energy deposited by the magnetic field. For this reason, and in order to keep the disk equations 
in strictly algebraically solvable form, we neglect the magnetic "frictional" heating. In all of these equations (except 
for the equation of state) we have neglected dimensionless coefficients of order unity. These factors have no impact on 
the form of the solutions for P, p, H, and T, and only a minor effect on the leading coefficients to the solutions. 

The four equations listed above, plus the equation for PH given in equation (|41|) can be solved algebraically in a 
manner analogous to that done for the original SS-disk equations to yield P, p, H, and T as functions of r, M, a, and 
/i. In solving these equations we have taken k to be given by Kramers opacity (re ~ 6 x 10 22 pT -3 ' 5 cm 2 gm _1 ) which 
is appropriate for most of the physical conditions found in our disk models, and the mass of the central star is fixed 
at 1.5 Mq. The solutions are: 

P~2x lO 5 ^ 9 / 1 ^ 72 Vo 21/8 ^ 17/2 ° dynes cm" 2 (45) 



3- — -2 




(42) 



H ~ 0.5 x lO^- 1 / 1 ^ 3 / 20 ^ 8 ^ 3 / 20 cm (46) 
T ~ 2 x lO^-^A^Vo 374 ^ 3710 K (47) 



p ~ 7 x lO-^- 7 / 10 ^'^ 1720 ^ 15 / 8 ^ 11 / 20 g cm- 3 (48) 

where M\q is the mass accretion rate in units of 10 16 gm sec -1 , and rio is the r adial distance in units of 10 10 cm. 
Analogous expressions have been derived by Matthews et al. (2005); see eq. (|B16[) in our Appendix B. Note that at 
r = r for any magnetic field value F(l, w) — /(?"o) = 0. Thus, our expressions for P, H, T, p (eqs. [45j 06J [47l [48] will 
vanish at r = r , just as in the Shakura-Sunyaev (1973) solution. This non-physical result is an artifact of the alpha 
viscosity prescription. Therefore, eqs. (|4"5l |4"51 |4"T1 14"8|) are not expected to be an accurate description of the disk close 
to r k ro. 

Eqs. |45 |) -([48 |) are derived under the assumptions that the accretion disks are optically thick, the radiative opacities 
are given by Kramers opacity, and gas pressure dominates over radiation pressure. The latter two assumptions can be 
changed to electron scattering opacity and/or a dominant radiation pressure, while still retaining analytic solutions. 
Some of these solutions, especially for values of M approaching Eddington, may be thermally unstable (e.g., the 
Lightman-Eardlcy instability; 1974), but this subject is beyond the scope of the present paper. 
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Fig. 9. — Surface density, S, vs. the radial distance (in units of r c ) for 20 values of the parameter £ in equal logarithmic steps 
over the range 0.1 < § < 10; with £ increasing from left to right. E, in units of g cm -3 , has been computed as the product of p 
and 2H, from eqs. I|48|l and (06j, respectively. The value of M used to evaluate S(r) has been arbitrarily fixed at 10 16 g s , but 
the curves scale simply as M 7 / 10 . The multiplicative constant "a" that multiplies £ is defined as r^ i ; this is done to allow the 
radial coordinate to be scaled in units of r c . 

In practice these equations could be utilized in the following way. For any accreting system which is thought to 
have a thin accretion disk, the accretion rate, and the mass, rotation rate, and magnetic moment of the central star 

2/7 

* ( T 1 (see eqs. [4], [5], [6]). This value of £ would 



h 2 /(Vgmm) 



would be combined to evaluate the parameter £ 
be used in eq. (fT7|) (or cq. |B7| . depending on the model for B^), to solve for ro, and for oj (eqs. [2D], [B10] ). With ro 
and lu in hand, eq. (|43|) (or eqs. [B18] and |B19j ) yields the function F(^ro/r,u>), which is required to solve for the 
pressure, density, and temperature in the disk via eqs. (I45l) -(l 



An illustrative example of disk radial profiles is shown as surface density plots in Fig. 9. Here we show £(r) for a 

but the curves 



range of values of the parameter £. The plots shown are for an arbitrary fixed value of M = 10 16 g s 1 



scale simply as M 7 / 10 . The multiplicative constant "s" that multiplies £ is defined as r c l0 , with r, 



c,10 



./(10 1 



this is done to allow the radial coordinate to be scaled in units of r c . In essence, up to a normalization factor, these 
are plots of y 3 ^ 2 F 7 ^ 10 (y,ui) versus w 2//3 y -2 , with each curve representing a different fixed value of u> (or, equivalently, 

x). 

8. SUMMARY AND CONCLUSIONS 

We have considered an analytic model of interaction between an accretion disk and a rotating magnetic field of the 
central star, in which the magnetic torques are distributed over a range of radii in the inner disk. Neglecting mass 
loss from the disk, we have solved the disk angular momentum equation to obtain a radial angular velocity profile 
smoothly matching the angular velocity of the star to that of the Keplerian disk. 

We have solved self-consistently for the inner radius of the viscous disk, and have found that it is within the corotation 
radius, ro < r c (§4). The magnetically dominated disk extends further in, with the maximum value of f2 attained at 
r\ < ro . We have found the viscous torques in the disk that follow from the specified model of magnetic torques (see 
§5). These, in turn, imply a luminosity profile of the disk, dL v - lsc /dr. While we have found the dependence of ro on 
the magnetic dipole for a specific model of magnetic field interaction with the disk, it is remarkable that most of the 
disk quantities presented in this paper need have no explicit dependence on the magnetic field, and depend directly 
only on M, ro, and r c . When expressed in a dimcnsionless form, i.e., in combinations such as r/r c or L/(GMM /r c ), 
the disk quantities depend on only one parameter, the ratio ro/r c , whose 3/2 power is the "fastness" u>. 

While most authors agree (for a fixed value of the magnetic dipole) as to the value of the mass accretion rate, 

M c ~ Mo = /i 2 r c 7//2 (GM) -1 / 2 , that corresponds to a termination radius of the disk that is more or less equal to the 
corotation radius, there is considerable theoretical disagreement as to what happens at lower accretion rates. A key 
question on which the authorities disagree is whether or not a strong magnetic field can prevent long-term accretion 
when the rate of flow of matter through the outer disk decreases below a critical value. 

It is possible that at very low mass-accretion rates, say M « Mo/10, the disk is expelled beyond corotation, and 
angular momentum is removed by outflowing matter (see the discussion in §1). Indeed, among the intermediate polars 
(in which there is clear evidence for the presence of both strong magnetic fields and an accretion disk) one system, 
AE Aqr, displays characteristics that may correspond to the propeller effect (Mikolajewski et al. 1996). Of course, for 
extremely strong stellar magnetic fields the disk may be disrupted completely, as in AM Her stars. 
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However, our results indicate that regardless of whether the disk is diamagnetic or not, and regardless of the detailed 
model of the disk-magnetosphere interaction, for mass accretion rates that are somewhat lower than the critical value, 
M < M c , there is a solution in which the disk terminates just within the corotation radius, r c . This is because the stellar 
torques must remove angular momentum from the disk to allow its termination, while outside r c the magnetosphcrc 
can only deposit angular momentum in the disk (Wang 1987). 

At these, and comparable, accretion rates (M ~ M c ), while removing angular momentum from the disk just inside 
the corotation radius, the magnetic torques may deposit angular momentum at r > r c , possibly at a substantial rate 
(depending on the torque model). This extra angular momentum must be transported outwards by viscous processes 
operating in the disk, leading necessarily to energy dissipation which increases the luminosity of the disk. Thus, in 
this spin-down regime, we report a qualitatively new result, which seems to have been overlooked in the literature. 
The angular momentum deposited in the disk by the pulsar torques leads to a substantial additional dissipation of 
energy in the disk (beyond the value expected from gravitational release alone) as angular momentum is transported 
outwards by viscous processes. This leads to a change in the total luminosity of the disk, L v ; s , as well as of its radial 
distribution (eqs. [26]. [28]). possibly quite dramatic at lower accretion rates. Consequently, previous estimates of the 
mass accretion rate of spinning-down pulsars may be wrong by up to a factor of several (Fig. 7), depending on the 
range of applicability of our model. 

The total contribution of magnetic origin to luminosity, already considered by Kenyon et al. (1996) in their study of 

T-Tauri stars is the magnetic torque power, L to t — GMM /(2ro), which includes the power L to t — Lvis released because 
of a mismatch in the rotation rate of the star and the parts of the accretion disk entrained by the magnetosphcrc 
(eqs. [3D], |33|). It is not clear a priori what fraction of the excess Ltot — L V j S is released inside the disk, and what 
fraction is released in non-thermal processes above the disk. 

We have quantified these processes within the context of two alternate specific models of distributed magnetic 
torques, both of which correspond to no shielding and no diamagnetism. Observations of torques as a function of 
luminosity may help in selecting the correct model of the magnetic field interaction with disk. For example, if the 
pulsed luminosity in X-ray pulsars is a good measure of the mass accretion rate, the observations discussed in Sj2] 
imply that in the spin-up phase the torque is proportional to M 13 , with j3 > 6/7. In our model I, at high accretion 
rates the torque is proportional to M 9 / 10 , while in the alternate prescription II (Appendix B) it is proportional to 
M 6 / 7 . However, the spin-up line for pulsars, eq. (|40|) . seems not to be very sensitive to the detailed model of magnetic 
interaction. The model of disk-magnctosphere interaction adopted in this paper does not quantitatively explain the 
reported rapid sign reversals of the torques in X-ray pulsars, with no apparent change in torque magnitude. However, 
already in our simple model the transition from spin-up to spin-down occurs over a fairly small luminosity range, 
and the spin-down torque quickly reaches an asymptotic value (Fig. 8), as the disk luminosity drops. This behavior 
has some qualitative resemblance to observations. It may be worthwhile to repeat the calculations for a partially 
diamagnetic disk. 

The results reported here include an analytic solution for the radial profile of angular velocity, both in the region 
where the viscous stresses do not vanish (which coincides with the region where the disk is Keplerian — this is the 
adopted ansatz), and inside the region where the viscous torques vanish. The main conclusions of the paper, such as 
that regarding the enhanced luminosity of the disk for spinning-down pulsars, are valid for any thin disk model (e.g., 
for any viscosity prescription). For a specific viscosity prescription, the a-disk, we have also computed the complete 
analytic radial profiles of the disk variables T, p, P, and H (§7). 

Note added in manuscript 

Since this paper was written the most extensive timing analysis on an accretion-powered millisecond X-ray pulsar 
(SAX J1808. 4-3658) has been reported by Hartman et al. (2007), who detect a long-term (i.e., over 7 years) spin down 
of the X-ray pulsar of v ~ —5.6 x 10~ 16 Hz s _1 . Perhaps of greater relevance to this paper, is the upper limit they 
set of \v\ < 2.5 x 10~ 14 Hz s _1 , during a typical X-ray outburst. We have used the X-ray outburst profile of 2002 
from Fig. 3 of Hartman et al. (2007) to numerically integrate our eq. (36) for the torque on the neutron star in order 
to determine the allowed range of neutron star magnetic moment, \i, that is consistent with the upper limit on the 
spin changes during a typical outburst. We find that fi< 3 x 10 26 G cm 3 , corresponding to a surface magnetic field 
strength of B < 3 x 10 8 G for a neutron star with a radius of 10 km. A similar value has been reported by Burderi 
et al. (2006). These are within the range of values cited by Hartman et al. (2007) for their constraint set during the 
outbursts; however, our model for the torques is substantially different than theirs. 

We thank Al Lcvinc, Ed Morgan, Ron Rcmillard, Jon Arons and Marek Abramowicz for helpful discussions. WK 
would like to acknowledge the hospitality of colleagues at the Kavli Institute at MIT. We are grateful to the anony- 
mous referee for drawing our attention to the work of Campbell and Matthews et al., as well as for other helpful 
comments. Antonia Savcheva's M.I.T. undergraduate senior thesis involved a number of studies related to this work. 
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APPENDIX 

APPENDIX A. PREVIOUS ESTIMATES OF THE MAGNETOSPHERIC RADIUS 

We briefly review a number of the earlier estimates of the disk termination radius that are found in the literature. 
Let us begin by asking at what value of a frozen-in poloidal magnetic field it would be energetically more favorable 
for plasma inside a flux tube to move out of the plane of the disk, rather than radially inwards. This would occur 
whenever the magnetic field is larger than the one given by the expression 

d2 

s = p > < A1 > 

where P is the plasma pressure in the inner disk. This is the same condition that was derived by Pringle & Rees 
(1972): "the radius down to which the gas can crush the stellar field (...) is the radius at which the effective pressure 
in the disk is equal to the external magnetic pressure" [assuming that this radius is within the corotation radius]. 
Neglecting screening of the field one substitutes the magnetic dipole value B = p/r 3 to obtain an estimate of this 
radius, tb through 

M P. (A2) 



87rr| 



To compare other workers' results, we recall that in a standard thin accretion disk all velocities scale as the Keplerian 
velocity vk = \/GAl/r 7 multiplied by a power of the dimcnsionless disk thickness H/r << 1 (Shakura & Sunyaev 
1973, Regev 1983, Kluzniak & Kita 2000). Thus, for example, the speed of sound scales as 

{P/p) 1 ' 2 ~c a ~(H/r)vK, (A3) 

while the radial drift velocity is 

v r ~ a(H/r) 2 v K . (A4) 

Here, p is the fluid density, and a is a dimensionless viscosity parameter. The radial angular momentum flux is then 

pv r VKT ~ apc 2 r ~ aPr. (A5) 

In the Ichimaru (1978) model, the disk is Keplerian everywhere except in a narrow boundary layer, in which all 
excess angular momentum is removed, so that at the base of the boundary layer, at r = ta, the fluid can come to rest 
in the stellar frame. The condition for hydrostatic equilibrium 

B 2 GMg /a n 

8^ = 

with cr ~ pi the surface density in the boundary layer of width I, translates into 

(A7) 

where Zbl ~ {H/r) 2 rA would be the thickness of a standard hydrodynamic boundary layer at this radius (Pringle 
1977). 

Scharlemann (1978) considers a completely diamagnctic disk and estimates the condition for stress balance between 
the disk and the screening toroidal current as 

B 2 ~ 21 P. (A8) 

We note that this is within the range of the other estimates discussed here, as 4tt < 21 < 87r. 

Arons (1993) considers a diamagnctic disk penetrated by the (enhanced) dipole over a region of radial extent 
comparable to the disk thickness, and writes the stress balance at the disk termination radius r = r m as 



B z B^ 

With v$ ~ vk, and B z ~ B^ this yields 



pVrV^. (A9) 



B 2 

— « aP. (A10) 

However, following the exact solution of Aly (1980) for the field configuration outside a diamagnetic disk, Arons takes 
the strength of the magnetic field to be enhanced by a factor of (r/H) 1 / 2 over the dipole value, yielding 

u 2 

- a(H/r m )P. (All) 
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Actually, Arons (1993) eliminates pv r with the help of the mass conservation equation for a steady disk, M 
to obtain a final condition on the inner radius of the disk in terms of the stellar mass accretion rate M: 



AirHrpVr, 
(A12) 

We have adopted this definition of convenient fiducial length. 

Wang (1987) obtained the same value (up to a factor of 2 1 / 7 ) by considering the condition that the magnetic stresses 
be able to remove enough angular momentum of the disk: 



r m = (GM)- 1/7 M- 2 /V /7 



B^B z r 2 = M-^(v K r) = ^Mv K . 



(A13) 



Wang (1987) expressly states that eq. (| A13|) holds only if the resulting termination radius is less than the corotation 
radius — otherwise, the central dipole would be unable to remove angular momentum from the disk. Again, with 
Btf, ~ B z , where B z is the unscreened and un-enhanced dipole value of the magnetic field, this would yield r = 2 1 / 7 r m 
as the disk termination radius. However, Wang (1987, 1995) uses a more complicated prescription for the value of B^, 
leading to an estimate very similar to our eq. (fl~7"|) in th e text. 

Clearly, as a - (H/r) << 1, the condition of eqs. (|A11[) and (|A12[) holds before the Pringle-Recs value of the magnetic 
field, i.e., eq. (|A1[) , can be attained. However, tb is smaller than r m only by a factor ~ (aH/r) 1 / 7 tn 0.4, where for 
the numerical estimate we took a as 0.01 and H/r as 0.1. 

In a very influential paper, Ghosh & Lamb (1979) suggested that a large region of the disk is threaded by a screened 
magnetic field of ~ 0.2 the strength of the external dipole, or less, and claim that the disk terminates at the radius 

r GL « 0.41 • 2- 1 / 7 r m . (A14) 

One concern may be that the value tql = 0.37r m is lower than the Pringle-Recs limit tb derived from eq. IjAip . 
and indeed, Ghosh & Lamb assume that matter leaves the disk vertically at a rate specified by an arbitrary "gating 
function." A number of authors (Pringlc & Recs 1972; Ichimaru 1978; Wang 1987; Spruit & Taam 1993; Rappaport, 
Fregau & Spruit 2004) restrict the termination radius to be within the corotation radius. 

In summary, the various authors differ by no more than a factor of 2 or 3 in their estimate of radius at which the 
central magnctosphcrc terminates a Kcplerian accretion disk, when this termination radius is within the corotation 
radius. Before the magnetic pressure becomes sufficiently strong to channel matter away from the symmetry plane 
of the disk, magnetic torques acting on the inner disk attain a value sufficient to remove angular momentum of the 
accreting matter at the rate corresponding to its advection at the inner edge of the disk, Mvkt. 

APPENDIX B. RESULTS FOR ALTERNATE PRESCRIPTION FOR B<j, 
Most of the results given in the text have been for the prescription where B$ is given by 



1 - — 



(Bl) 



for all radii, i.e., r < r c and r > r c . Here we give the principal analytic results for the alternate prescription for Bl 
that we have considered: 



BV~Bi for r>r 



B% 1 ~ B^us/n) ~ -B z ( 1 - ^ ) for r < r c 



(B2) 
(B3) 



and for r < r c the last term in eq. (TTS")) is then replaced by an alternate form of the magnetic torque 

Blr 



1 - -4r 



4ttF v n 

For the region r > r c the prescription remains the same as for B\. The corresponding vertically integrated r 
component of the viscous stress-energy tensor, analogous to eq. (I24|) in the text is: 



(B4) 




for r < r c 



for r > r c 



The boundary condition for the inner edge of the viscous accretion disk, analogous to eq. (TT7| . is (Wang 1995): 



1 = 2 



ro 



7/2 , 




(B5) 
(B6) 

(B7) 
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The asymptotic forms of eq. (|B7[) are 

r /r c ~2 2 / 7 £ for £«1 , 
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(B8) 



r /r c ~ 1 - 7/2 for £ » 1 



(B9) 



which is the same as eq. (fTTj)) . For a system with fixed /U, uj s , and M, eq. (|B8|) implies that r$ cx r m cx M 2 / 7 in the 
high mass accretion rate limit of eq. (|B7|) . The counterpart of eq. (|2"0j) is 



(r m /ro) 7/2 = [2(l-w)] -1 , 



(BIO) 



or 2Mo(l — u>) = Mu) 7 ' 3 , and the same limits obtain as for eq. (j20|) : for M — * 0, u — > 1; for fj, — * 0, u> — > 0. However, 
now r m /r — ► 1/2 2 / 7 (and not 0) as /i — ► 0. The expression for dL v - ls /dr corresponding to eq. l|2"6")) in the text is given 
by: 



dL vis _ 3GMM 
dr 2r 2 

dL vis _ 3GMM 
dr 



i (^y^ 2 



1+1 



rV 
roJ 



for r < r c 



3 , r .3/2 / r 2 X 3 / 2 

2[ — ) +2( — ) +6( 



r r c 



2r 2 

for r > r c . 

The total torque on the central star, analogous to eq. (|36l) in the text is given by 

,2 



\ 3 / x 3/2 

7V. / V r r . 



For R s — > and M/0, this reduces to 

r cs = MyjGMr [7/6 - (4/3)w + (l/9)w 2 ] /(l - w) 
(Wang, 1995). The dimcnsionlcss torque corresponding to eq. (|37|) is given by gu(oj), with 



9ii{u) 



1 



1 3w-2w 2 
21 1 - w 



(Bll) 



(B12) 



(B13) 



(B14) 



(B15) 



The first root is at w = [(12- v / 102)/2] ps 0.9502 w (0.9666) 3/2 , quite close to the root of gi in prescription I (eq. [59]): 
20/21. 

Finally, the factor F, corresponding to eq. (T4"2"|) is 



F= 1 



7/2 



7/2 



2 ' ' 3(1/ r .3 



ro 



for r < r c 



3/2 



''(] 



- 3 



for r > r c 



(B16) 



(B17) 



To within a factor of two, the magnetic torque of eq. (|B4[) was adopted by Matthews et al. (2005) for all r, accordingly, 
to within the same factor of two, our eq. (|B7[) is the same as their eq. (26), although its interpretation is different (see 
our discussion of eq. [TT] in § I4.2[) . After the sub stitut ion of parameter values f3 — > /1 2 / (71V GM), 7 = 7/2, eqs. (30-33) 
of Matthews et al. (2005) coincide with our eq. (|B16j) in the region of overlap of t he to rque prescription in their paper 
and this Appendix (r < r c ). However, for r > r c , we obtain a different result, eq. (|B17[) . because Matthews et al. have 
adopted a different prescription for this regime. 

After using eq. (|B10j) we obtain the function F(y/ro/r,ui) that appears in the a-disk equations 

F(y,w) = [6(l-w)]~ 1 (6- 7y + y 7 - 6w + 8wy - 2wj/ 4 ) , for u 1/3 < y < 1 , (B18) 
F(y,cu) = [18(1 - lo)]- 1 (18 - 2ly - 3y 7 - 18w + 2Alu V - 2uj 2 y + 2(J- 1 y w ) , 

for < y < uj 1/z . (B19) 

Again F(l,u>) = 0, and the zero magnetic field limits of eqs. (|B18[) . (|B19[) are F(y,0) = 1 — y — (y — y 7 ) /6, and 
F(y, 0) = 1 — y — (y + t/ 7 ) /6 — ► 1, respectively. Note that in view of the domain restriction y < cj 1 / 3 , the last term in 
eq. (jB19[) . i.e., w _1 y 10 , is not divergent when uj — > 0. 
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APPENDIX C. A MAGNETICALLY DOMINATED KEPLERIAN ACCRETION DISK 

The specific prescription for that we are using, is incompatible with a magnetically dominated Keplerian disk, i.e., 
the solution to eq. (|16p is not Keplerian. This prescription is equivalent to assuming the following form for magnetic 
diffusivity: rji ~ rHQ. Campbell (1992) uses a different form, T)c = vaH, where the Alfven speed va ~ B$/ '(inp) 1 / 2 . 
His work cannot be directly compared with ours, because we restrict ourselves to a discussion of a thin accretion disk. 
Eq. (53) of hydrostatic equilibrium in Campbell (1992) neglects gravity, but for a gravitationally stratified disk with 
sound speed, ~ (P/p) 1 ^ 2 ~ QrH, equal to the Alfven speed va, Campbell's diffusivity would differ from ours by a 
factor of H/r: r/c — T)\(H/r). [In fact, in the Campbell (1992) solution, H ~ r, so that in terms of their numerical 
values rjc ~ 771 .] 

Inspired by Campbell's approach, in this Appendix we discuss the solution for a gravitationally stratified, Keplerian 
accretion disk that is dominated by magnetic torques. We take B z = pr~ 3 , M = iirrv r pH, r/c = HrH 2 , 

S /(4 7 rp) 1 / 2 = v A = n K H, (CI) 



B^ = rHB z (fl K - ui s )/r, c = B z {r/H){\ - w./Hk), (C2) 



and consider the region within the corotation radius, fix > w s . 
The angular momentum conservation equation 



implies that 



The height of the disk now follows directly from eq. (|C2[) . while the density from eq. (|C1 



M± (fSl K ) = ^ (C3) 
or r 



Mn K r 2 1/2 

^ = — 2p] — r ^ ' 



Mtl K r 4 V tt K ) V fix/ 



Note the divergence in density at the corotation radius. 
The radial velocity starts from zero at the corotation radius, and rapidly increases inwards: 

-1 „ ,,-15/2 ' w - 



v r oc (rffp)- 1 oc r~ lb ' 2 1 - ^ . (C7) 

Since all of the angular momentum is removed locally from the disk by the magnetic torques, the torque on the star 
from an annulus at r is (1/2)Q,k Mr dr, see eq. (|C3p . and the power applied by the annulus to spin up the star is 
(l/2)MQjfW s rdr. At the corotation radius all the power of the disk is used to torque the star, and so the disk is 
dark; for smaller radii the disk becomes brighter as the relative amount of angular momentum removed decreases. The 
radiated flux is then 

**■<- 5* (•-&)■ < C8 > 

the first term being the rate of energy change of the fluid over unit area 

2irrdr \ 2 r J 4tt k ' V ' 

the second term is the spin- up power drained from the disk per unit area. This expression, eq. (|C8j) . agrees with the 
Joule heating formula derived by Campbell (1992): 

2F 7 = 2r 2 B 2 z (n K - u s ) 2 H/ric. (CIO) 



The total torque on the star is 



Mtt K rdr = Ml K (r c ), (Cll) 



assuming Rg -C r c . Here, Ik = ^K r2 is the Keplerian specific angular momentum. 

We conclude that with a suitable prescription for a magnetic field, it is possible to construct a magnetically dominated 
accretion disk that formally has Keplerian orbital velocities. In fact, however, the trajectories of the fluid would not 
be close to Keplerian orbits for all radii — in addition to its azimuthal motion, the fluid suffers rapid radial inflow. 



Magnetically Torqucd Accretion Disks 



19 



Because the radial velocity i s co ntrolled directly by the prescribed magnetic field — at steady accr etion rate, eq. (|C7|) 
ultimately follows from cq. (jC2[) and from the equation of vertical hydrostatic equilibrium (eq. |Clj ) — the solution 
cannot be considered valid for all r < r c until its consistency with the radial equation of motion has been checked. 

In the main body of this paper, we find a solution for the angular velocity profile in the magnetically dominated 
region which is not Keplerian at all (Fig. [2j. The counterpart of eq. (|C8j) for the region r < tq is 



T = 



M 
2nr 



(n 



n 2 



— (nr 2 ) (n - u s ) 



(C12) 



This is the power released per unit area of the disk over and above the mechanical power used to torque up the star. 
The first term in p arentheses in eq. (|C12|) is the power used to accelerate the fluid radially. For ft = Qk, eq. i|C12[l 
reduces to cq. (IC8j) . 
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